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A strategy is proposed to design the broadband gain-doped epsilon-near-zero (GENZ) 
metamaterial. Based on the Milton representation of effective permittivity, the strat- 



egy starts in a dimensionless spectral space, where the effective permittivity of GENZ 
metamaterial is simply determined by a pole-zero structure corresponding to the op- 
erating frequency range. The physical structure of GENZ metamaterial is retrieved 
from the pole-zero structure via a tractable inverse problem. The strategy is of 
great advantage in practical applications and also theoretically reveals the cancel- 
lation mechanism dominating the broadband near-zero permittivity phenomenon in 
the spectral space. 
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Materials with near-zero permittivity (ENZ) have been vastly explored due to their 
anomalous electromagnetic properties and unique applications, such as directive antenna 1 , 
subwavelength imaging system 2 ' 3 , optical nanocircuit 4 , radiation phase pattern converter 5 , 
reflectionless waveguide 6 ' 7 , field energy confinement 8 , and electromagnetic transparency and 
cloaking design 9 ' 10 . In nature, such materials are already available as noble metals and polar 
dielectrics. Characterized by the Drude or Drude-Lorenz model 11 ' 12 , the real part of their 
permittivity approaches to zero near the plasma frequency. However, the broadband ENZ 
material is required regarding to the practical applications. To manufacture the broadband 
ENZ materal, the basic idea is to make the metal-dielectric composite according to the ef- 
fective medium theory 13 . But such mixture possesses a high loss in optical frequency range 
due to the Kramers-Kronig relations 14 , which limits the application scenario. To overcome 
such disadvantage, gain sources are suggested to reduce the loss 15 , and due to the strong 
local field enhancement in ENZ material, the gain sources can even provide large effective 
gain than when used alone 16,17 . 

Based on former study about the graded metal-dielectric mater ials 18 ~ 20 , we have promoted 
a strategy to design an anisotropic broadband ENZ metamaterial by applying the Bergman 
representation of effective permittivity 21,22 in our previous work 23 . However, a high loss 
exists in the operating frequency range due to the absence of gain sources. In addition, the 
previous strategy requires analytical expressions of permittivities of different components 
to determine the physical structure of the broadband ENZ metamaterial. That limits its 
practical applications, where permittivities are usually numerically given 24 . Therefore, we 
promote another strategy to design the anisotropic broadband gain-doped epsilon-near-zero 
(GENZ) metamaterial base on Milton representation of effective permittivity 25-27 in this 
work. 

The strategy is carried out in a dimensionless spectral space, in which the effective per- 
mittivity of broadband GENZ metamaterial is characterized by a series of pole- zero pairs, 
initially arranged under lossless condition with respect to a properly determined operating 
frequency range. Then with proper variations on the initial zeros (zero- variation, for short), 
the real part of effective permittivity can be near- zero in the operating frequency range when 
loss is taken into account. Finally, via a tractable inverse problem, the physical structure 
of broadband GENZ metamaterial can be achieved from the varied zeros and poles. In 
contrast to the previous one, this strategy can be carried out numerically with respect to 
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the experimental data, and due to the gain sources, the loss can be effectively reduced. In 
the following, the strategy will be described in detail, and it is worth noting that in order to 
simplify the description, permittivities are still characterized by analytical formulae, which 
is not necessary. 

Depicted in Fig. 1, the broadband GENZ metamaterial is a flat iV-layer stack. Each 
layer is composed by two components: metal-gain capsules doped in a dielectric host. The 
metal-gain capsule is a metallic core covered by a gain shell. The permittivity of metallic 
core follows the Drude model 

e m (u) = e 00 1 (1) 

CO [co + 17) 

where is a high frequency permittivity determined to match the experimental data in 
visible region, and 7 is the damping factor. The gain shell is composed by two kinds of gain 
molecules doped in dielectric host, with an effective permittivity described by the Lorentz 
model 15 

e g (u) = e r + — — — + -= — — , (2) 

where e r is the permittivity of dielectric host, (k%, k 2 ) are the density factors with respect to 
the number density of each gain molecule, and (woi, W02) are the center emission frequencies 
of each gain molecules with the emission frequency line-widthes (uj, a 2 ). For convenience, 
frequencies (u, u i, U02), damping factor 7, and line-widthes (cri a 2 ) are all normalized by 
the plasma frequency u p of metallic core. Regarding the simple mixing rule, the effective 
permittivity of metal-gain capsule can be described as 

e 1 (u)=pe g (u) + (l—p)e m (u) 1 (3) 

where p is the volume fraction of gain shell. 

For simplification, the permittivity of dielectric host of each layer is set to be 

e 2 = e r . (4) 

Besides, the thickness of ith-layer is denoted as d\ and normalized by the total thickness of 
GENZ metamaterial. Therefore, the following summation rule 

N 



E^ 1 ( 5 ) 



is held. 

For simplification, the effective permittivity of ith layer is described by the simple mixing 
rule, as 

e®{u,) = f i e 1 {u) + 0—f t )e 2 , (6) 

where is the volume fraction of metal-gain capsules in the ith layer. By introducing a 
s-parameter 

Eq. (6) can be rewritten as a function of s-parameter 

ef{s)= e 2 (l- £). (8) 
Therefore, the effective permittivity of GENZ metamaterial is 

N 



£e[s = e 2 



. -i 

di 



(9) 



t 1 - f*/ s 

According to this equation, a series of proper volume fractions (fi) and thicknesses (di) can 
make the effective permittivity of GENZ metamaterial, in the normal direction of the stack, 
be near- zero in a wide frequency range. 

On the other hand, the effective permittivity of GENZ metamaterial can also be described 
by the Milton representation as 

N 

ee(s)=e 2 H S -^, (10) 

1=1 

where the pole- zero (si-zi) pairs are restricted as 

^ Si < Z\ < s 2 < z 2 < ■ ■ ■ < s N < z N < 1. (11) 

Therefore, the effective permittivity of GENZ metamaterial is simply determined by a series 
of pole-zero pairs in contrast to Eq. (9). 

Based on the foregoing text, the design of the broadband GENZ metamaterial can be 
simply arranged as follows: 

1. Determine a proper operating frequency range, in which not only the real part of effec- 
tive permittivity of GENZ metamaterial can be near- zero, but also the loss (imaginary 
part of effective permittivity of GENZ metamaterial) can be effectively reduced by 
gains. 



2. Determine a series of pole-zero pairs in order to make the real part of effective per- 
mittivity of GENZ metamaterial be near-zero in the operating frequency range. 



3. Determine the volume fractions and thicknesses with respect to the pole-zero pairs via 
an inverse problem, in order to obtain the physical structure of GENZ metamaterial. 

In the following, these steps are introduced in detail, and an example is depicted in Fig. 2. 

In order to reduce the loss, the operating frequency range is dominated by the effective 
permittivity of gain shell and the s-parameter. The imaginary part of effective permittivity 
of gain shell dominates the effect of gain molecules. Therefore, it is proper to locate the 
operating frequency range between the emission frequencies (w i, W02) to make a good use 
of the gain molecules. However, the dramatic variations of the effective permittivity of gain 
shell at the emission frequencies should be avoid. Therefore, the density factors (ki, K2) 
should be carefully chosen to make the effective permittivity varies smoothly between the 
emission frequencies, just as depicted in Fig. 2(a). 

On the other hand, the imaginary part of s-parameter should also be considered. Just 
as indicated in Fig. 2(b), it is clear that between the emission frequencies, the imaginary 
part of s-parameter has two zero points (u)&, oj C 2). Therefore, the final operating frequency 
range, denoted as (u a , Ub), should be localized between the two zero points. 

The second step is to determine a series of pole-zero pairs with respect to the operating 
frequency range, and it can be finished as follows: 

(1). Determine the initial zeros and the corresponding poles. 

i). This step starts with the determination of zeros. According to Eqs. (7), (10), and 
(11), it is clear that all zeros should be real functions of frequency. Therefore, it 
is simple to define the zeros as the real part of s-parameter 




(12) 



where the frequency Ui is determined as follows 



u>i = u a + (u b - cu a )(i - 1)/(N - 1) (for 2 = 1,2,..., N), 



(13) 



and named as zero-frequency. 
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ii). Regarding Eq. (11), the corresponding poles can be determined as follows: 

(for i = l), 

(14) 

Zi-x + Zi)/2 (for i = 2, 3,..., N). 

It is worth noting that due to the simple mixing rule, i.e., Eq. (8), the first pole 
should be zero. 

Based on the definition of zeros and poles, the effective permittivity exactly equals 
zero at all zeros and extends to infinity at all poles under lossless condition, which is 
depicted in Fig. 2(c) in orange. However the result changes when loss is taken into 
account, i.e., the dashed curves in Fig. 2(c). Therefore, the initial zeros should be 
revised under lossy condition. 

(2). Zero- variation under lossy condition. 

To vary the initial zeros under lossy condition, the following equation is considered 



»M*(w))U = » 



' N i \ i 
£2n fM^i 



(for i — 1,2, ... , N), (15) 



where the varied zeros are denoted as z[, and the values of frequency uji and poles 
Si are from Eqs. (13) and (14) respectively. Note that it is difficult to solve this 
equation directly. However, based on the values of initial zeros, i.e., Eq. (12), it is can 
be numerically solved. After the variation, the real part of effective permittivity will 
equals to zero at the varied zeros, which is indicated as solid curve in Fig. 2(c). 

According to Eqs. (5), (9), and (10), the inverse problem reads 

E^v=nF7- 1 ' (16) 
i=i s h i=\ * 

with respect to the varied zeros z\ and the poles Sj. The right hand side of Eq. (16) can be 
expanded into 

i=i 1 i=i 1 

Therefore, the volume fraction of metal-gain capsules in ith layer is just /j = z\, and the 
thickness of ith layer is = QiJ z[. After all these steps, the physical structure of GENZ 



metamaterial is finally determined. The detailed derivation of Eqs. (16) and (17) is given in 
the Appendix. 

Figure 2(c,d) displays a schematic example about the design. Figure 2(c) displays the 
variations of effective permittivity of GENZ metamaterial with respect to frequency in dif- 
ferent conditions. In lossless condition, dt(e e (oj)) exactly equals zero at every initial zero- 
frequency and extends to infinity at the frequency corresponding to the poles. When loss is 
taken into account, 9R(e e (u;)) becomes continuous in the operating frequency range. However, 
without zero-variation, it is not near-zero in the operating frequency range. In addition, a 
finite 9f(e e (a;)) is also introduced due to the loss. In contrast, with the varied zeros, 9ft(e e (u;)) 
becomes near-zero in the operating frequency range. Yet the varied zeros also effect $s(e e (uj)), 
due to the Kramers-Kronig relations. 

Figure 2(d) offers a comparison between the results with and without gain. Clearly, 
^s(e e (u)) can be dramatically reduced by gains. In this example the average reduction on 
the loss can be about 50%. 

Finally, it is worth noting that by increasing the number of layer, reducing the width of 
operating frequency range, or properly increasing the volume fraction of gain shell in the 
metal-gain capsule, better result can be obtained. 

The zero- variation strategy is not only effective in designing GENZ metamterial, but also 
reveals the physical mechanism of broadband ENZ phenomenon, i.e., the cancellation mech- 
anism in electric displacement. Moreover, due to the Milton representation, the cancellation 
mechanism is vividly indicated in the spectral space. 

Figure 3 takes a schematic example to explain the cancellation mechanism in spectral 
space. Consider a 2-layer stack, the effective permittivity of each layer follows Eq. (8), and 
the volume fraction follows j\ > fa for convenience. Initially, there is only one pole-zero pair 
in the spectral space for the first stack only, i.e., sx = and z\ — fa, according to Eqs. (8) 
and (10). By adding the second layer, it amounts to add an additional pole-zero pair (s 2 and 
z-i) in the spectral space. On the left hand side of the additional pole S2, the polarization 
is out of phase with the response of the pole S\, thus yielding a zero Z2 between the poles. 
Therefore, the out-phase polarization can cause the cancellation of electric displacement, 
leading to a near- zero effective permittivity. By properly arrange a series of pole-zero pairs, 
the broadband ENZ can be achieved. Contrary to the physical space, the spectral space has 
no concern with physical structure, thus it offers a more clear picture about the cancellation 
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mechanism. 

In this work, a strategy based on the Milton representation of effective permittivity is 
proposed to design the broadband GENZ metamaterial. The application of Milton repre- 
sentation not only simplifies the design procedure, but also reveals the physical mechanism 
dominating the broadband ENZ phenomenon. Meanwhile, the potential applications of this 
strategy in practical design is also implied, since it is can be carried out numerically. 

Furthermore, not only limited in the design of flat multi-layer broadband GENZ meta- 
material, the strategy also offers a fundament to design broadband GENZ metamaterial 
in other geometrical structures, and material with other permittivity in a wide frequency 
range. 



APPENDIX 

The derivation of Eqs. (16) starts with Eq. (9), which can be written as 

N a 

Regarding Eq. (5), the above equation can be rewritten as 

1=1 

On the other hand, Eq. (10) gives 



fi/ s ) ~^ s ~fi 



N 

e i=l 1 

Finally, with respect to the varied zeros z' i: the inverse problem, Eq. (16), reads, 

N fH N «« «, 

EJi a i _ TT b ~ b i 1 

Regarding the derivation of Eq. (17), it can be numerically carried out by applying Mathe- 
matica 8 of Wolfram Research, or other similar computational softwares. 
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FIGURE CAPTIONS 



FIG. 1. (Color on line) The broadband GENZ metamaterial is a flat iV-layer structure 
(a). Each layer is composed by metal-gain capsules (b) doped in dielectric host. The metal- 
gain capsule is a metallic core covered by a gain shell. The permittivity of metallic core 
is denoted as e m (oS). The gain shell is a shell of dielectric host doped with two kinds of 
gain molecules, and its effective permittivity is denoted as e g (uj). In addition, the effective 
permittivity of metal-gain capsule is denoted as Si(u), calculated by the simple mixing rule, 
and the permittivity of dielectric host is set as Ei = e r . The thickness of ith layer is denoted 
as di, restricted by the summation rule. 

FIG. 2. (Color on line) Figure 2 indicates a schematic example about the GENZ meta- 
material design with the following parameters (i) high frequency permittivity = 1.0, (ii) 
damping factor 7 = 0.02, (iii) permittivity of dielectric host 82 = e r = 1.0, (iv) emission 
frequency u i = 0.2 and Uq\ = 0.3, (v) emission frequency line- width o\ = 0.02 and a 2 = 0.1, 
(vi) density factor K\ = 0.02 and k 2 = 0.015, and (vii) volume fraction of gain shell p = 0.6. 
Figure 2 (a, b) displays the variation of effective permittivity of gain shell and the imaginary 
part of s-parameter, which determine the location of operating frequency range. Figure 2(c) 
denotes the effective permittivity of a 10-layer GENZ metamaterial with respect to the op- 
erating frequency range 0.22 ~ 0.27 determined by Fig. 2(a, b). Under lossless condition, 
3^(£ e ( w )) I lossless exactly equals to zero at all zero-frequencies, but extends to infinity at the 
frequencies relating to the poles. With respect to the initial zeros, 5R(^ e (u;)) | initial is not 
near-zero in the operating frequency range, but increases slightly with respect to frequency. 
After the zero-variation, 9ft(£e(w))| varic d presents a good result in the operating frequency 
rang. Furthermore, it is clear that with respect to the near-zero result of 3fJ(e e (o;))| varicd , 
Q^eeO^-O) I varied ^ s a little greater than 9 : (£ e (a;))| initial , due to the Kramers- Kronig relations. 
Figure 2(d) gives a compare between the results with and without gains. Clearly, the gain 
does not affect the real part of effective permittivity too much in the operating frequency 
range. However, the loss is dramatically reduced by the gains and the average reduction is 
about 50% in this example. 

FIG. 3. (Color on line) Figure 3 gives a schematic explanation about the cacellation 
mechanism in spectral space, which causes the broadband ENZ phenomenon. It indicates 
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the effective permittivity of a 2-layer stack (red-solid), and the permittivity of each layer 
(red- and orange-dashed). Consider only the first layer, there is only a pole-zero pair (si and 
Zi). By adding the second layer, it amounts to add an additional pole-zero pair (s 2 and z 2 ) 
in the spectral space, thus yielding the out-phase polarizations by the sides of the additional 
pole s 2 . The out-phase polarizations can cause the cancellation of electric displacement, 
leading to the ENZ phenomenon. Therefore, with a series properly arranged pole- zero pairs, 
the ENZ phenomenon can occur in a wide frequency range. 
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FIG. 2. L.Sun and K.W.Yu 
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FIG. 3. L.Sun and K.W.Yu 
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